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The  paper  concerns  solution  manifolds  of  nonlinear  parameter-dependent 


equations  (1)  F(u,A)  *  y^  involving  a  Fredholm  operator  F  between  (infinite- 

"  T3  /**.u*- 

dimensional)  Banach  spaces  X  s  Z  x  ^  and  Y,  and  a  finite-dimensional 

jt.rr'  Li-*-* 

parameter  space  y.  Differential -geometric  ideas  are  used  to  discuss  the 
connection  between  augmented  equations  and  certain  one-dimensional  submanifolds 
produced  by  numerical  path-tracing  procedures.  Then,  for  arbitrary  (finite) 

s.  /**  0  t  ^ 

dimension  of  j(,  estimates  of  the  error  between  the  solution  manifold  of  (1) 

■  and  its  discretizations  are  developed.  These  estimates  are  shown  to  be 
applicable  to  rather  general  nonlinear  boundary-value  problems  for  partial 
differential  equations. 


1.  Introduction 


Considerable  attention  has  been  directed  toward  nonlinear  parameter- 
dependent  equations  of  the  form 

0.1)  F(u,X)=y0, 

where  F  is  a  mapping  from  a  Banach  space  X  =  Z  x  A  to  a  Banach  space  Y, 

Z  represents  a  state  space,  and  A  an  m-dimensional  parameter  space, 

1  <  m  <  oo.  Under  appropriate  conditions  on  F,  the  set  of  regular  solutions 
of  (1.1)  forms  an  m-dimensional  manifold  in  X  and  it  is  of  interest  to  obtain 
both  analytical  and  computational  information  about  this  manifold. 

Since  Z  and  hence  X  are  usually  infinite-dimensional,  any  computational 
examination  of  (1.1)  requires  the  introduction  of  finite-dimensional  approxi¬ 
mations.  Such  approximations  raise  questions  about  the  resulting  error  between 
the  solution  manifolds  of  (1.1)  and  its  discretized  versions. 

The  estimation  of  this  error  between  solution  manifolds  has  seen  relatively 
limited  activity,  and  most  of  that  has  been  directed  toward  the  case  of  a  one¬ 
dimensional  parameter  space  A  [3],  [4],  [6].  Since  all  numerical  procedures 
for  analyzing  the  solution  manifold  of  an  equation  of  the  form  (1.1)  involve 
tracing  one-dimensional  submanifolds,  consideration  of  this  special  case  is, 
in  a  sense,  reasonable.  However,  it  is  generally  not  easy  to  construct  a  correct 
picture  about  an  m-dimensional  manifold  from  information  obtained  on  certain  one¬ 
dimensional  submanifolds.  Accordingly,  it  is  certainly  desirable  to  derive  error 
estimates  for  the  general  case  of  an  m-dimensional  parameter  space. 

This  is  the  topic  of  the  present  paper.  The  definition  of  the  errors  under 
consideration  depends  critically  on  the  choice  of  the  local  coordinates  on  the 
manifolds.  After  some  preliminaries  in  Section  2,  these  local  coordinates  are 


the  topic  of  Section  3.  Then  in  Sections  4,5,6  we  consider  in  some  detail 
the  idea  behind  the  numerical  procedures  used  to  trace  paths  on  a  manifold. 
Included  here  are  the  role  played  by  augmented  equations  and  their  use  in 
characterizing  points  on  the  manifold. 

Then  in  Section  7  we  extend  the  error  estimates  developed  in  [4]  to 
the  general  situation  of  m-dimensional  solution  manifolds.  These  results 
allow  for  the  study  of  approximations  of  such  manifolds  directly,  without 
recourse  to  paths.  The  arguments  of  Section  7  depend  on  the  validity  of  a 
certain  stability  condition  whose  implications  and  various  reformulations  are 
the  topic  of  Section  8.  Finally,  two  examples,  in  Section  9,  illustrate  the 
theory  and  show,  in  particular,  that  it  applies  to  a  very  wide  class  of  boundary 
value  problems  for  partial  differential  equations. 
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2.  Prel iminaries 

Throughout  this  paper,  we  shall  use  the  following  assumption: 

F:  ECX  +  y  is  a  Cr-Fredholm  mapping  with  r  _>  1  and  index 
(A)  m  1  from  an  open  subset  E  of  a  Banach  space  X  into  a 

Banach  space  Y. 

As  usual,  a  point  x  £  E  is  a  regular  point  of  F  if  DF(x)  is  surjective. 
For  any  yQ  e  F(E),  our  interest  centers  on  the  regular-solution  set 

(2.1)  M  =  {x  e  E:  x  regular,  F(x)  *  yQ), 
of  the  equation 

(2.2)  F(x)  *  yQ. 

The  fundamental  result  concerning  the  structure  of  M  is  contained  in  the 
following  theorem  proved  in  [4]. 

Theorem  2.1 :  For  any  yQ  e  F(E) ,  the  regular  solution  set  M  is  a  relatively 
open,  m-dimensional  Cr-manifold  in  X. 

The  notion  of  parametrization  introduced  in  [4]  extends  easily  to  the 
manifold  M.  However,  we  modify  the  terminology  somewhat  in  order  to  be 
consistent  with  the  customary  language  of  differentiable  manifolds  and  bi¬ 
furcation  theory. 

For  any  xQ  e  M,  let  Tx  M  denote  the  tangent  space  of  M  at  xQ.  In 

o 

our  setting,  we  can  identify  T  M  with  ker  DF(x  ).  Let  X  =  V©T  be  any 

*o  0 

splitting  of  X  such  that  dim  T  =  m  and  VAT  M  =  (0).  With  any  such 

*o 


splitting,  it  is  always  possible  to  define  an  isomorphism  A:  Y  V  of 
Y  onto  V.  In  particular,  the  splitting  X  =  V©T  M  with  any  comple- 

o  1 

mentary  subspace  V  of  T  M  and  the  isomorphism  A  *  (DF(x  ) | V)  has 

xo  0 

important  applications  in  the  reduced-basis  technique  (see  [5]). 

The  following  theorem  is  an  immediate  generalization  of  the  corre¬ 
sponding  result  in  [4].  Its  proof  involves  an  application  of  the  implicit 
|  function  theorem. 

Theorem  2.2:  Let  yQ  e  F(E),  X  *  V  ©T  be  a  splitting  as  above  at  a  given 
point  xQ  e  M,  and  A:  Y  -*■  V  an  isomorphism  of  Y  onto  V.  Then  there 
exist  an  open  ball  JC  T  with  0  e  J,  an  open  neighborhood  UCX  of  xQ, 
and  a  unique  Cr- function  n:  J  -*■  Y  such  that  n(0)  *  0  and 

M  n  U  =  (x  e  X:  x  *  x(t)  =  xQ  +  t  +  An(t),  t  e  J}. 

This  theorem  justifies  the  use  of  the  word  parametrization  in  [4], 
However,  since  we  wish  to  reserve  the  word  parameter  for  another  concept,  we 
have  to  carry  the  result  of  Theorem  2.2  somewhat  further. 


Corollary  2.3:  Under  the  conditions  of  Theorem  2.2,  there  exist  an  open 
ball  JQCLJ  with  0  e  Jfl  and  an  open  neighborhood  UQ  c  U  of  xQ  such 
that  t  x(t)  -  x0  +  t  +  An(t)  Isa  Cr-diffeomorphism  of  JQ  onto  M  r\UQ 


I 


Proof:  From  F(x(t))  =  0  for  t  e  J,  it  follows  that  DF(xQ)x'(0)  *  0  and 


hence,  since  T  and  V  are  complementary  subspaces,  that  x'(0):  T  -*•  T  M 


Is  an  isomorphism.  Let  W  be  any  complementary  subspace  of  T  M  in  X 

xo 


and  define 


<f>:  J  +  W  CX  -*■  X,  <p(  x)  *  x(t)  +  w,  x*t+weJ  +  W. 


Then  we  have 


4>'(0)C=x'(0)t+co,  C  =  t+  u)cT  M  ®  W  , 

xo 

and  thus  ^'(O)  is  an  isomorphism  of  X.  By  the  inverse  function  theorem, 

<p  maps  a  neighborhood  of  0  Cr-diffeomorphically  onto  a  neighborhood  of 
0.  Hence,  because  of  <j>(t  +  0)  =  x(t),  the  result  follows. 

The  inverse  of  the  mapping  x:  JQ  HflU0  is  a  chart  or  coordinate 

mapping  at  xQ  of  the  manifold  M  and  such  mappings  are  often  said  to 

define  a  system  of  local  coordinates  on  M.  For  that  reason,  we  call  a 

splitting  X  *  V  ®  T  such  that  dim  T  =  m  and  VHT  M  =  {0}  a  coordinate 

xo 

splitting  for  M  at  xQ  e  M,  and  we  refer  to  T  as  a  coordi nate  space. 


3.  Parameters  as  Coordinates 


It  often  happens  in  applications  that  certain  quantities  are  naturally 
identified  as  parameters.  This  means  that  there  is  an  intrinsic  splitting 
X  =  Z©A  of  X,  where  A  is  an  m-dimensional  parameter-space  and  Z 
represents  a  state  space.  Such  a  splitting  of  X  will  be  called  a  parameter 
splitting.  It  is  natural  to  attempt  to  use  the  parameter  space  A  as  the 
coordinate  space  T  of  a  coordinate  splitting,  and  the  question  arises  when 
this  is  possible. 

In  order  to  provide  an  answer  to  this  question,  we  suppose  for  the  re¬ 
mainder  of  this  section  that  a  parameter  splitting  X  =  Z®  A  is  available. 
Moreover,  for  a  given  yQ  e  F(E)  let  M  denote  the  corresponding  regular- 
solution  manifold  (2.1)  and  xQ  any  point  of  M. 

The  suitability  of  A  as  the  coordinate  space  t  of  a  coordinate 
splitting  depends  on  the  subspace 

(3.1)  Z  =  ZAT  M. 

o 

If  ZQ  =  {0},  then  clearly  X  =  Z©A  is  a  coordinate  splitting  for  the  mani 

fold  M  at  x  .  But  if  Z  +  {0}  then  T  M  has  a  nontrivial  component 
oo  xQ 

in  Z  and  hence  we  cannot  use  all  of  A  as  the  coordinate  space  T. 

Minimally  we  should  exchange  ZQ  for  a  part  of  A  not  meeting  Tx  M. 

o 

To  this  end,  let  n  be  the  natural  projection  of  X  onto  A  along  Z 
and  define  the  space 

(3.2)  A  =  nT  M. 

*n 

Moreover,  the  following  nota*.»>r  wi  be  convenient:  If  W  is  any  Banach 
space  and  WQCW  a  closed  subspace  which  splits  W,  then  W  ©  WQ  shall 
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denote  any  closed  complementary  subspace  of  WQ  in  W. 

Some  important  properties  about  the  subspaces  ZQ  and  Aq  are  con¬ 
tained  in  the  following  lemma. 

Lemma  3.1 : 

(i)  dim  A„  *  m  -  dim  Z  . 

o  o 

(ii)  Aq  =  Aft  (Z  +  Tx  M). 

0 

(ill)  z©  A  =  Z  +  T  M. 

0  Xo  | 

(iv)  [(Z0Zo)©(A0Ao)]nIx  M-  {0}. 

0 

Proof:  The  proof  involves  only  elementary  notions  about  subspaces.  Clearly 
(i)  needs  no  elaboration.  To  prove  (ii),  let  A  e  A  .  Then  A  =  z  +  u  for 
Sc..a  u  e  T  M  and  any  z  e  Z,  whence  A  e  AO  (Z  +  T  M).  Conversely,  if 

A  A 

0  0 

the  latter  inclusion  holds  for  some  A  then  A  =  z  +  u  for  some  z  e  Z, 

u  e  T  M  and  A  =  IIA  *  JIu  e  A  .  In  (iii)  the  containment  Z  ©A  C  Z  +  T  M 
*o  0  0  xo 

is  obvious  in  view  of  (ii).  To  show  the  reverse,  let  x  =  z  +  u,  z  e  Z, 

u  e  T  M.  Then  there  exist  z1  e  Z,  A  e  A  such  that  x  =  z'  +  A,  whence 

xo 

A  =  (z-z1)  +  u  e  Z  +  T  M  and  thus  A  e  A  by  (ii)  or  x  e  Z©A  .  Finally, 

xQ  o  o 

let  u  belong  to  the  set  on  the  left  of  (iv)  and  u  =  z  +  A,  z  e  Z,  A  e  A. 

Then  A  =  -z  +  ue  Z  +  TM  =  Z©A.  whence  A  e  An  and  A  e  A„0  ( A  ©  A  )  =  {0}. 

Xq  0  0  0  0 

As  a  resul  t  we  have  u  =  zeZ©ZCZ  and  therefore  u  e  l  O  T  M  =  Z 

o  xo  0 

which  implies  that  x  e  ZQ0  (Z©  ZQ)  *  {0}. 

I 

The  above  facts  can  be  combined  to  give  the  desired  theorem  on  coordinate 


spl ittings. 


Theorem  3.2:  In  addition  to  assumption  (A),  let  X  =  Z©  A  be  a  param¬ 
eter  splitting  of  X,  yQ  e  F(E),  and  xQ  a  point  of  the  corresponding 

regular-solution  manifold  M.  Let  Z  =  Z(T  T  M  and  A  =  nT  H,  where 

xo  0  xo 

II  is  the  projection  onto  A  along  Z.  Then  X  =  [(Z©  ZQ)  ©  (A©  Aq)] 
©  [Zq  ©  Aq]  is  a  coordinate  splitting  for  M  at  xQ. 

This  theorem  is  illustrated  in  the  next  section. 
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4.  Parameters  and  Integral  Manifolds 

All  currently  available  numerical  methods  for  analyzing  a  solution  mani¬ 
fold  of  a  parametrized  equation  involve  the  computational  trace  of  one¬ 
dimensional  submanifolds  by  means  of  a  continuation  process.  These  one¬ 
dimensional  submanifolds  are  usually  defined  by  combinations  of  the  natural 
parameters  with  one  degree  of  freedom  and  the  specification  of  a  starting 
point.  In  differential -geometric  terms,  this  means  that  we  compute  specific 
integral  manifolds  of  certain  1 -distributions  on  the  manifold. 

Before  we  define  these  notions  explicitly,  an  example  to  illustrate 
the  situation  may  be  in  order.  Consider  the  buckling  of  a  spring  system 
discussed  in  [  9  ;  pp.  301-305],  After  a  suitable  scaling,  the  total  energy 
of  the  system  is 

2  1  2 

U(p,q,A,v,y)  =  (1-p)  +  ^  Y(2q)  +  2Apcosq  +  vpsinq, 

where  p,q  are  displacements  (the  state  variables)  and  the  parameters  A,v,y 
represent  applied  forces  and  a  spring  constant,  respectively.  The  spring 
constant  y  is  intrinsically  positive.  As  a  result,  the  equilibrium  equation 
is 


(4.1) 


F(x) 


-2(1 -p)  +  2Acosq  +  vsinq 


=  0, 


4yq  -  2Apsinq + vpcosq 


where  F:  EC3R  ->-lR  is  a  C  -mapping  on  E  *  {x  =  (p,q,A,v,y)  elR  :  y  >  0}. 
It  is  easly  verified  that  all  points  of  E  are  regular  for  F  and  that 
0  e  F(E).  Consequently,  the  equilibria  of  the  system  correspond  to  points  on 
the  three-dimensional  manifold  M  =  (x  e  E:  F(x)  =  0}.  The  natural  parameter 


splitting  X  =  Z©A  associated  with  the  parameters  A,v,y  is  obtained  by 
setti ng 


(4.2)  Z  =  {(p»q»0,0,0)T) ,  A  =  {(0,0,A,v,y)T}. 


As  indicated  before,  we  introduce  now  some  combination  of  parameter 
value'  with  one  degree  of  freedom.  The  simplest  approach  is  to  fix  two 
of  the  above  three  parameters,  say  v  =  0,  Y  =  This  is  equivalent  to 
introduction  of  the  augmented  equation 


(4.3) 


G(x)  = 


5  4  oo 

where  G:  EC1R  -*-]R  is  still  of  class  C  on  E,  but  now  not  all  points 

of  E  are  regular  for  G.  For  example,  at  xq  =  (^- ,  0,  j  ,  0,  ^)T  e  E, 

DG(xq)  has  rank  3  whereas,  of  course,  the  rank  of  DF(xQ)  is  2. 

The  point  xQ  =  (j,0,p0,g-)T  is  an  example  of  a  point  where  the 
natural  parameter  splitting  (4.2)  of  the  original  problem  is  not  a  coordinate 
splitting  and  Theorem  3.2  comes  into  play.  For  this  particular  xQ,  we  have 


and  hence  that 


(4.4) 


T  M  =  ker  DF(x  )  =  {(p,q,-p,0,y)  }, 
o 


Z  =  ZnT„  M  =  ((0,q, 0,0,0)  }, 
0  *o 


A_  =  ITT  M  =  {(0,0, M,y)'}. 
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By  Theorem  3.2,  we  see  that 

X  =  [(Z©Z0)©(A©A0)]©[Z0®  aq] 

(4.5) 

=  {(p,0,0,v,0)T}  ©  {(0,q,X,0,y)T} 

is  one  possible  coordinate  splitting  for  M  at  xQ. 

Generally,  let  M  denote  any  finite-dimensional  Hausdorff  manifold  of 
class  Cr,  r  1 .  We  again  denote  the  tangent  space  of  M  at  a  point  x  e  M 
by  TxM  and  write  TM  for  the  tangent  bundle.  Recall  that  a  vector  field 
of  class  Cp,  1  <  p  <  r,  on  a  relatively  open  subset  MQ  C  M  is  a  Cp-mapping 
Mq  TM  such  that  £(x)  e  TxH  for  each  x  e  MQ.  Moreover,  a  1 -distribution 
A  of  class  Cp  on  Mq  is  defined  by  the  following  two  properties: 

(i)  A  is  a  mapping  A:  Mq  -*  TM  such  that  Ax  is  a  one-dimensional 
subspace  of  T  M  for  each  x  e  M  . 

A  U 

(ii)  For  each  xQ  e  Mq,  there  exist  an  open  neighborhood  U  of  xQ 

in  X  and  a  vector  field  M0  O  U  TM  of  class  Cp  such  that 
A  *  span  U(x)}  for  all  x  e  M  A  U. 

A  U 

A  one-dimensional  submanifold  N  of  MQ  is  an  integral  manifold  of  A  on 
Mq  if  for  every  x  e  N  the  tangent  space  TXN  is  equal  to  A  .  For  an 
introduction  to  these  concepts  see,  for  example,  [11]. 

For  our  simple  example,  let  EQCE  be  the  open  set  of  all  regular 
points  of  G  and  set  MQ  =  MnEQ.  Then  A:  MQ  -*■  TM,  Ax  =  ker  DG(x), 
x  e  Mq,  is  a  mapping  which  associates  with  each  regular  point  x  e  MQ  a 
one-dimensional  subspace  of  the  tangent  space  TxM.  In  this  case,  it  happens 
to  be  possible  to  define  a  global  vector  field  on  MQ  which  satisfies  the 
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second  condition  of  a  1 -distribution.  In  fact,  for  any  x  e  EQ  a  vector 

c 

£(x)  eR  is  uniquely  defined  by 


(4.6) 


DG(x)£(x)  =  0, 


I|5(x)M2  -  l. 


det 


'DG(x)' 

.S(*)T. 


>  0, 


M  1 

and  it  is  readily  verified  that  the  mapping  x  ►+•  £(x)  is  of  class  C 

on  E„  and  hence  a  C  -vector  field  on  fi.  This  shows  that  A  is  indeed 
o  o 

p.l 

a  1 -distribution  of  class  C  on  M^. 

o 

The  point  xQ  *  (j* 0*  \ » 0,  j)T  turns  out  to  be  a  singular  point  of  the 
vector  field  £.  In  fact,  for  the  value  of  y  *  g-,  the  bifurcation  set  in  the 
Av-plane  consists  of  a  butterfly  and  touching  dual  cusps  (see  Figure  1  and 
[9]).  The  dual -cusp  point  corresponds  to  xQ. 

The  global  definition  (4.6)  of  the  vector  field  is  possible  only  in 
the  finite-dimensional  case.  However,  we  do  not  need  a  determinant  to  define 
an  orientation  of  a  local  vector  field  in  a  neighborhood  of  each  x  e  Mq. 

In  other  words,  it  may  be  expected  that  augmented  equations  of  the  form  (4.3) 
always  define  1-distributions  on  certain  relatively  open  submanifolds  MQ  of 
M.  The  numerical  methods  mentioned  earlier  are  then  designed  to  compute 
specific  integral  manifolds  of  such  1-distributions.  The  connection  between 
augmented  equations  and  certain  1 -distributions  is  taken  up  in  detail  in  the 
next  section.  Once  augmented  equations  are  associated  with  1-distributions, 
they  can  then  be  used  to  classify  an  arbitrary  point  of  M  and  to  provide 
an  alternate  formulation  of  Theorem  3.2.  This  is  the  topic  of  Section  6. 

Our  example  shows  that  the  singular  behavior  at  the  point  xQ  *  (j,  0,  j ,  0,  gOT 
arises  only  because  we  considered  a  particular  parameter  combination  and  the 
corresponding  1 -distribution.  Otherwise,  the  point  is  regular  for  F.  More¬ 
over,  Instead  of  fixing  v  ■  0,  y  *  ^  and  letting  A  vary  as  we  did  in  (4.3), 


Ei 


t' 

► ' 


K 


we  may  fix  v  =  y  *  g-  and  allow  v  to  vary.  Then  we  are  led  to  the 
augmented  equation 

F(x) 

(4.7)  G(x)  =  X 

Y 

For  this  new  G,  DG(xq)  has  rank  4  so  that  xQ  now  becomes  a  regular 
point  of  G.  This  opens  up  the  possibility  of  considering  numerical  methods 
which  avoid  the  singularity  entirely  by  working  with  other  more  suitable 
1 -distributions  such  as  that  defined  by  the  augmented  equation  (4.7). 


5.  1- Distributions  and  Augmented  Equations 

In  the  last  section,  we  discussed  the  idea  behind  continuation  methods 
for  analyzing  solution  manifolds  of  parametrized  equations  and  how  parameter 
combinations  lead  to  1-distributions  and  one-dimensional  integral  manifolds. 

In  this  section,  we  consider  the  relationship  between  such  integral  manifolds 
and  certain  augmented  equations. 

Throughout  this  and  the  next  section,  let  M  denote  the  m-dimensional 

M 

C  -manifold  (2.1)  of  the  equation  (2.2)  and  assume  that  X  =  Z©A  is  a 
parameter  splitting  of  X.  A  combination  of  the  natural  parameters  with  one 
degree  of  freedom  corresponds  to  the  choice  of  a  one-dimensional  subspace 
SCA  which  defines  the  remaining  degree  of  freedom.  We  call  S  a  reduced- 
parameter  space  and  any  spl i tti ng  X  =  W  ©  S  with  Z  C W  a  reduced- parameter 
splitting.  For  example,  in  (4.3)  the  reduced -parameter  space  is 
S  =  {(0,0,A,  0,0)T),  whereas  in  (4.7)  S  =  {(0,0,0,v,0)Th 

r— 1 

Suppose  now  that  A:  MQ  -*•  TM  is  a  1 -distribution  of  class  C  on 
some  open  set  M  c  M  for  which  IIA„  =  S  for  all  x  e  M,.  As  in  Section  3, 
n  denotes  the  projection  of  X  onto  A  along  Z.  Clearly,  the  condition 
IIAX  s  S  is  equivalent  to  AXCZ©S.  We  call  a  1 -distribution  with  this 
property  a  1 -distribution  with  respect  to 

1-distributions  with  respect  to  S  can  be  discussed  conveniently  in 
the  context  of  certain  augmented  forms  of  the  equation  (2.2).  For  any  given 
reduced-parameter  space  S,  there  exist  linear  operators  L:  A  with 

ker  L  =*  S.  With  any  choice  of  such  L,  we  define  the  augmented  mapping 

(5.1)  G:  ECX  +  Y  x  ]Rm-1 ,  G(x)  =  (F(x),LHx),  x  e  E, 


and  with  it  the  augmented  equation 
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(5.2)  G(x)  =  (y0,Lna0), 

where  aQ  e  X  is  as  yet  arbitrary.  For  example,  in  (4.3)  and  (4.7) 
L:  A  1R^  is  defined  by 


o 

o 

o 

o 

1 _ 

o 

o 

o 

o 

II 

and  L  = 

1 

o 

o 

o 

o 

1 _ 

0  0  0  0  1 

respectively. 

A  number  of  useful  facts  about  (5.1)  are  collected  together  in  the 
following  technical  lemma. 

Lemma  5.1:  Let  x_  e  M.  Then: 

-  o 

(i)  06(xq)  is  a  Fredholm  operator  of  index  l. 

(ii)  DF(xq)|Z,  and  hence  DG(xQ)|W,  is  a  Fredholm  operator  of  index  0. 

(iii )  ker  DG(xq)  =  Tx  MH  (Z®  S). 

o 

(iv)  Z  =  ZDT  M  =  wn  ker  DG(x  ). 

U  Aq  o 

(v)  S  O  ITT  M  »  {0}  If  and  only  if  1  =  ker  DG(x  ). 

*0  0  0 

Proof:  Note  that 

DG(xq)  =  (DF(x0),Ln)  =  (DF(xo),0)  +  (O.Ln), 

where  (DF(xQ),0):  X  Y  x  is  a  Fredholm  operator  with  index  m  -  (m-1 )  *  1 

and  (0 , LIT) :  X  -►  Y  x  ]Rm'^  is  a  compact  operator.  This  implies  that  (i)  holds 
(see  [10;  p.  114]).  The  injection  j:  Z  -*•  X,  jz  *  z,  z  e  Z,  is  a  Fredholm 
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operator  of  index  0  -  m  =  -m.  Therefore,  DF(xQ)|Z  =  DF(xQ)j  is  a 

Fredholm  operator  with  index  m  -  m  =  0  ([10;  p.  Ill])  which  proves  (ii). 

Now  (iii)  follows  directly  from  the  definition  of  G  and  the  fact  that 

u  e  ker  DG(x  )  if  and  only  if  u  e  ker  DF(x  )  =  T  M  and  nu  e  ker  L  =  S. 
o  0  XQ 

Furthermore,  (iv)  is  a  consequence  of  (iii)  and  Z  C  W.  Finally,  to  obtain 
(v),  we  note  that 

so  a  =  {oi <=>  n(T  m  0(ze  S))  =  {0} 

0  xo 

<=>  nker  DG(xq)  =  {0}  by  (iii) 

<=>  ker  DG(xQ)C  Z 

<=>  ker  DG(>^)  =  lQ  by  (iv). 


This  completes  the  proof. 

Note  that  none  of  the  assertions  of  the  lemma  depends  upon  the  particular 
choice  of  L  in  the  definition  of  G.  A  consequence  of  part  (iii)  is  the 
following  connection  between  1 -distributions  with  respect  to  S  and  augmented 
equations. 

Theorem  5.2:  Let  X  =  W©S  be  a  reduced-parameter  splitting  and  A:  MQ  -*•  TM 
a  1-distribution  with  respect  to  S  of  class  C  on  a  relatively  open  set 
H  CM,  Moreover,  let  G  be  an  augmented  function  (5.1).  Then  A  C.  ker  DG(x) 
for  any  x  e  MQ  and  hence  any  integral  manifold  N  of  A  on  MQ  is  a 
solution  manifold  of  the  augmented  problem  (5.2),  where  aQ  is  a  given  point 


of  N. 


So  far,  we  began  with  a  1 -distribution  with  respect  to  S  and  then 
related  it  to  an  augmented  equation.  In  practice,  the  starting  point  is 
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usually  an  augmented  equation  which  is  then  used  to  characterize  and 
determine  a  1-distribution  with  respect  to  S.  This  was  the  procedure  we 
followed  in  the  example  of  the  previous  section.  How  this  can  be  accomplished 
in  general  is  our  next  goal. 

As  before,  let  X  *  W©S  be  a  reduced-parameter  splitting  of  X  and 

G  an  augmented  function  (5.1).  We  denote  by  EqC  E  the  open  set  of  all 

regular  points  of  G  and  set  MQ  =  M  A  Eq. 

Theorem  5.3:  The  mapping  A:  Mq  -*■  TM,  Ax  =  ker  DG(x),  x  e  Mq,  is  a  1- 

distribution  with  respect  to  S  of  class  C  on  Mq.  Moreover,  if  r  >_  2 > 

then  for  any  given  aQ  e  MQ  the  regular-solution  set  of  the  augmented 
equation  (5.2)  is  an  integral  manifold  of  A  on  MQ. 

Proof:  By  Lemma  5.1 (i),  DG(x)  is,  for  any  x  e  Mq,  a  Fredholm  operator  of 
index  1  and  since  DG(x)  is  surjective  we  must  have  dim  DG(x)  =  1. 
Consequently,  the  mapping  A  is  well-defined  and  by  Lemma  5.1 (iii)  we  have 

Ax  =  TxM  O  (Z©S)  for  x  e  MQ.  Therefore,  it  remains  to  verify  only  the 

second  property  of  a  1 -distribution.  By  the  theory  of  Fredholm  operators  [2], 

[10],  for  xQ  e  M  ,  there  exist  a  vector  uQ  e  X  and  a  continuous  linear 

functional  <pQ  e  X*  such  that  Ax  =  ker  DG(xQ)  *  span{uQ},  <{>0(u0)  *  1, 
and  X  =  ker  <J>Q  ©  Ax  .  The  mapping 

H:  E  x  X  +  Y  x  ]R,  H(x,u)  =  (DG(x)u,l-4>Q(u)),  (x,u)  t  E  *  X, 

is  of  class  C  and  satisfies  H(x0*u0)  =  0.  As  in  the  proof  of  part  (i) 
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of  Lemma  5.1,  the  partial  derivative 

DuH(xo’uo)*  =  (DG<x0)-‘°)  +  (0.-*0O) 

is  the  sum  of  a  Fredholm  operator  of  index  0  and  a  compact  operator  and 
hence  is  itself  a  Fredholm  operator  of  index  0.  If  DuH(xq,u0)v  =  o, 
then  v  e  ker  4>qAAx  =  {0};  that  is,  DuH(xq,u0)  is  an  isomorphism.  By 
the  implicit  function  theorem,  there  exist  a  nieghborhood  U  of  xQ  in  X 

r_l 

and  a  C  -mapping  £:  U  -*•  X  such  that 

H(x,C(x))  =  (DG(x)£(x).l-<Do(£{x)))  =  (0,0),  x  e  U. 

Therefore,  £:  MQn  U  TM  is  a  local  vector  field  and  Ax  =  span{£(x)} 

for  x  e  Mq A  U.  This  completes  the  proof  that  A  is  a  1 -distribution  with 

1 

respect  to  S  of  class  C  .  The  second  part  follows  from  the  standard 
existence  and  uniqueness  theory  for  flows  since  we  assumed  r  >_  2  (for 
example,  see  [  8]). 

Theorem  5.3  guarantees  that  the  regular-solution  set  of  an  augmented 
equation  can  be  used  to  determine  an  integral  manifold.  But  an  augmented 
function  (5.1)  is  defined  on  all  of  E;  that  is,  an  augmented  problem  can 
also  be  considered  at  nonregular  points  of  the  augmented  function.  We  consider 
this  topic  in  the  next  section. 


20 


6.  Characterization  of  Points  Using  Augmented  Functions;  Coordinate  Splittings 


We  continue  the  development  begun  in  the  previous  section  and  consider 
how  an  augmented  function  can  be  used  to  characterize  the  nature  of  an 
arbitrary  point  on  the  regular-solution  manifold  M  of  equation  (2.2)  with 
respect  to  a  given  parameter  splitting  X  =  Z©A  and  reduced-parameter 
space  S.  Theorem  5.3  applies  only  on  the  submanifold  Mq  of  M  consisting 
of  the  regular  points  of  the  augmented  function.  Nevertheless,  any  augmented 
function  is  defined  everywhere  on  M  and,  as  noted  earlier,  that  observation 
enables  us  to  use  an  augmented  function  to  characterize  different  types  of 
points  on  M. 

A  general  examination  of  the  possible  types  of  points  on  M  is  not  our 
intent.  Instead,  we  consider  three  commonly  occurring  cases  which  arise 
frequently  in  practice.  The  following  definition  is  based  upon  fairly 
standard  terminology. 


Definition  6.1 ;  Let  X  *  W©S  be  a  reduced-parameter  splitting  and  G  any 
corresponding  augmented  function  (5.1).  A  point  xq  e  M  is  a 

(i)  nonsingular  point  of  G  if  dim  ker  DG(xQ)  =  1  and 
wn  ker  DG(xQ)  =  {0}; 

(ii)  limit  point  of  G  if  dim  ker  DG(xQ)  =  1  and 
ker  DG(xQ) C  W;  and 

(iii)  simple  critical  point  of  G  if  dim  l,.er  DG(xQ)  =  2  and 
dim  (WO  ker  DG(xQ)}  =  1. 


Note  that  this  definition  depends  only  on  W  and  S  and  is  otherwise  in- 


dependent  of  the  choice  of  G.  In  fact,  the  definition  involves  only  W 
and  ker  DG(xq)  and,  in  view  of  Lemma  5.1 (iii),  ker  DG(xQ)  depends  only 
on  S  and  not  on  the  particular  choice  of  the  linear  operator  L:  A  -"lRm"^. 
Also,  we  remark  that  nonsingular  points  and  limit  points  of  G  are  regular 
points  of  G. 

For  the  three  cases,  we  may  characterize  completely  the  nature  of  a 

point  x  e  M  in  terms  of  the  subspaces  Z  =  Z  fj  T  M,  A  =  nT  M,  and 
o  0  xo  0  xo 

S.  In  the  following  discussion,  we  assume  the  setting  of  Theorem  3.2,  (5.1), 
and  Definition  6.1.  We  also  assume  that  r  >  2  in  assumption  (A)  (in  order 
to  apply  Theorem  5.3). 

We  begin  with  the  case  1Q  =  {0}.  Then  from  Lemma  5.1 (ii)  and  (iv), 
it  follows  that  xQ  is  a  nonsingular  point  of  G.  By  Theorem  5.3  (with 
aQ  =  xQ)  and  Theorem  3.2  (applied  to  G),  X  -  W®  S  is  a  coordinate 
splitting  for  the  regular-solution  submanifold  of  the  augmented  problem 

(6.1)  G(x)  =  (y0»LnxQ) 

at  xQ.  Moreover,  by  Lemma  5.1 (iv)  and  Theorem  3.2  (applied  to  F), 

X  =  Z®  A  is  a  coordinate  splitting  for  the  manifold  M  at  xQ. 

Next  we  assume  that  dim  ZQ  =  1  and  SC\  AQ  =  {0}.  From  Lemma  5.1  (v), 
it  follows  that  xQ  is  a  limit  point  of  G,  and  Theorems  5.3  and  3.2  imply 
that  X  =  [(W©  ZQ)  ©  S]  ©  ZQ  is  a  coordinate  splitting  for  the  regular- 
solution  submanifold  of  the  augmented  problem  (6.1)  at  xQ.  (Note  that  the 
use  of  ZQ  =  ker  DG ( xQ )  as  the  coordinate  space  is  the  basis  for  the 
Lyapunov-Schmidt  procedure  [12].)  To  extend  this  coordinate  splitting  for 
the  submanifold  to  a  coordinate  splitting  for  the  full  manifold  M,  we 
combine  Lemma  3.1  (i)  with  Lemma  5.1  (v)  to  conclude  that  S©  A  =  A  and 


then  use  Theorem  3.2.  This  shows  that  X  =  [(Z  ©  Z)  ©  S]  ©  [Z„©  A  ]  is 
a  coordinate  splitting  for  M  at  xQ. 

Finally,  let  dim  ZQ  *  1  and  SO  Aq  /  {0}  (that  is,  S  0  AQ  =  S). 

Again  from  Lemma  5.1,  it  follows  that  xQ  is  a  simple  critical  point  of  G. 
Although  DG(xq)  remains  a  Fredholm  operator  of  index  1,  its  null  space  is 
now  two-dimensional  and  xQ  is  not  a  regular  point  of  G.  Nevertheless,  xQ 
is  still  a  regular  point  of  F  and  the  manifold  M  does  have  coordinate 
splittings  at  xq.  In  particular,  by  Theorem  3.2  once  more, 

X  =  [(Z©  Z  )  ©  (A©  A  )]  ©  [Z  ©  A  ]  is  a  coordinate  splitting  for  M  at 
x  .  It  is  important  to  observe  that  now  S  is  not  an  acceptable  choice  for 
the  complementary  subspace  A©  A  . 

We  summarize  the  above  results  cn  coordinate  splittings  in  our  next 
theorem. 

Theorem  6.2:  Assume  the  conditions  and  notation  of  Theorem  3.2,  (5.1),  and 
Definition  6.1. 

(i)  If  ZQ  =  {0},  then  xQ  is  a  nonsingular  point  of  the  augmented 

function  G;  X  =  W©S  is  a  coordinate  splitting  for  the  regular- 
solution  submanifold  of  the  augmented  problem  (6.1)  at  xq,  and 
X  =  Z©A  is  a  coordinate  splitting  for  the  regular-sol ution  mani¬ 
fold  M  of  (2.2)  at  xQ. 

(ii)  If  dim  ZQ  =  1  and  SO  AQ  =  {0},  then  xQ  is  a  limit  point  of 

the  augmented  function  G;  X  =  [(W  ©  ZQ)  ©  S]  ©  lQ  is  a  coordinate 
splitting  for  the  regular-solution  submanifold  of  the  augmented 
problem  (6.1 )  at  xQ,  and  X  =  [(Z  ©  ZQ)  ©  S]  ©  [ZQ  ©  AQ]  is  a  co¬ 
ordinate  splitting  for  the  regular-solution  manifold  M  of  (2.2)  at 
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(iii)  If  dim  ZQ  =  1  and  SO  A0  /  {0},  then  xq  is  a  simple  critical 

point  of  the  augmented  function  G;  X  =  [(Z©  ZQ)  ©  (A©  AQ)]  ©  [Z  ©  AQ] 
is  a  coordinate  splitting  for  the  regular-solution  manifold  M 
of  (2.2)  at  xQ.  Here,  S  is  not  an  acceptable  choice  for  A©Aq. 

Thus,  from  a  knowledge  of  the  subspaces  ZQ,  Aq,  and  S  alone,  we  are 
able  to  characterize  the  nature  of  the  point  xQ  with  respect  to  a  combination 
of  parameters.  In  Theorem  6.2(iii),  it  is  not  possible  to  identify  the  type 
of  simple  critical  point  because  that  would  require  an  examination  of  second 
and  higher  derivatives.  However,  it  is  interesting  to  note  that  the  coordinate 
splitting  for  the  manifold  M  is  independent  of  the  type  of  simple  critical 
point  at  xQ;  in  other  words,  it  is  not  necessary  to  know  the  type  of  simple 
critical  point  when  choosing  a  coordinate  splitting  for  M. 

To  illustrate  this  theorem,  we  return  to  the  problem  of  the  buckling  of 
a  spring  system  mentioned  in  Section  4.  With  xQ  =  (^-,  0,  ,  0,  g-)T,  we 

calculated  ZQ  and  Aq  in  (4.4).  Therefore,  according  to  Theorem  6.2,  since 
dim  ZQ  =  1 ,  xq  is  either  a  limit  point  or  a  simple  critical  point  of  any 
augmented  function  G.  The  augmented  function  G  of  (4.3)  is  associated  with 
the  reduced-parameter  space  S  =  {(0,0,X,0,0)T)  for  which  S  fl  AQ  f  (0).  By 
Theorem  6.2(111),  xQ  must  be  a  simple  critical  point  of  that  augmented 
function  and  (4.5)  is  a  coordinate  splitting  for  the  regular-solution  manifold 
of  (4.1)  at  xQ.  On  the  other  hand,  the  augmented  function  G  of  (4.7)  is 
associated  with  the  reduced-parameter  space  S  =  {(0,0,0,v,0)T)  for  which 
S  A  Aq  *  (0).  By  Theorem  6 . 2 ( i i ) ,  xQ  is  a  limit  point  of  that  augmented 
function;  X  =  {(p,0,X,v,y)T}  ©  { (0 ,q ,0 ,0 ,0 ) ^}  is  a  coordinate  splitting  for 
the  regular-solution  submanifold  of  the  augmented  problem  (4.7)  at  xQ,  and 
(4.5)  again  serves  as  a  coordinate  splitting  for  the  regular-solution  manifold 
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7.  Finite- Dimensional  Approximations  and  Error  Estimates 

We  consider  now  the  fundamental  problem  of  determining  suitable  approxi¬ 
mations  to  the  regular-solution  manifold  M  of  the  equation  (2.2)  and  of 
obtaining  corresponding  estimates  for  the  error.  As  it  turns  out,  the  dis¬ 
cussion  in  [4  ]  applies  essentially  unchanged  to  the  case  of  an  m-dimensional 
manifold.  Therefore,  we  merely  introduce  the  notation  and  state  the  principal 
results  in  our  current,  more  general  setting  and  refer  to  [4  ]  for  the  proofs 
which  are  easily  extended  to  this  case. 

As  always  we  assume  that  the  information  in  (A)  is  given.  To  formulate 
an  appropriate  approximate  problem  for  (2.2),  we  remark  that  in  applications 
it  is  the  state  variables  and  not  the  natural  parameters  which  are  discretized. 
With  that  in  mind,  suppose  that  a  parameter  splitting 

(7.1 )  X  =  Z  ©  A  (dim  A  =  m) 

of  X  is  available  and  that  Z  and  Y  are  related  through  an  operator  Q 
as  follows: 


(7.2) 


Q  e  L(X,Y) ,  ker  Q  =  A, 

Q | Z  is  an  isomorphism  of  Z  onto  Y. 


A  discretization  of  the  problem  (2.2)  is  specified  by  a  collection  (P^:  h  >  0} 
of  finite-rank  projections  Ph  e  L(Y)  converging  strongly  to  the  identity  on 
Y;  that  is. 


lim  P.y  =  y,  y  e  Y. 
h-0  n 


(7.3) 
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The  projections  and  the  operator  Q  determine  finite-dimensional  sub¬ 
spaces 

(7.4)  Yh  =  PhY*  zh  =  (Qlz)"lYh.  xh  =  Zh®A‘  T 

and  out  interest  is  directed  towards  the  regular- sol ution  manifold 

(7.5)  ^  =  {x  e  Eh:  x  regular,  Fh(x)  =  yQh> 
of  the  equation 

$4 

<7-6>  Vx>  '  *oh-  yoh  *  Vo- 

where 

(7.7)  Fh;  EhC  xh  -  Yh,  Fh(x)  -  PhF(x),  xeEh  =  Enxh. 

We  refer  to  the  information  contained  in  (7.3)-(7.7)  as  a  basic  discretization  F3 

of  the  problem  (2.2). 

It  is  easier  to  compare  a  discretized  problem  (7.6)  with  the  original 
problem  (2.2)  if  both  are  formulated  on  the  same  spaces.  A  convenient  way 
of  doing  this  is  to  extend  the  discretized  mapping  F^  to  all  of  EC  X  by 
defining 

(7.8)  Fh:ECX-Y,  Ph(x)  =  (I-Ph)Qx  +  Ph(F(x)-yQ) ,  x  e  E, 

where  I  denotes  the  identity  on  Y. 

Some  important  properties  of  and  its  connection  to  a  discretized 
problem  (7.6)  are  summarized  in  the  next  lemma. 


i  re 

a  A 


Lenina  7.1 :  Fh  is  a  Cr-mapping  with  the  following  properties: 

(i)  Fh(x)  =0  for  x  e  E  if  and  only  if  x  e  Eh  and  Fh(x)  =  yQh. 

(ii)  DFh(x)XhcYh,  x  e  E. 

(iii )  ker  DFh(x)  c^.  x  e  E. 

( i v )  DFh(x)  c  L(X,Y)  is  a  Fredholm  operator  of  index  m  for  x  e  E. 

(v)  PhDF(x)Xh  *  Yh  for  some  x  e  E  implies  that  x  is  a  regular  point 
of  P. 

To  compare  the  regular-solution  manifolds  M  of  (2.1)  and  of  (7.5), 
let  xQ  be  any  point  of  M.  Moreover,  let  X  -  V&T  be  a  coordinate 
splitting  for  M  at  xQ  and  A:  Y  V  an  isomorphism  of  Y  onto  V.  The 
relationship  between  the  coordinate  splitting  and  the  above  basic  discretization 
is  expressed  in  the  following  stability  condition: 

(s)  HDPh(x0)Ay|  |  1  6I  lyl  I*  y  e  Y»  h>0  sufficiently  small , 

where  S  is  a  positive  constant  independent  of  y  and  h.  A  detailed  dis¬ 
cussion  of  the  implications  of  this  stability  condition,  including  equivalent 
formulations,  is  deferred  to  the  next  section. 

The  stability  condition  (S)  combined  with  generalized  versions  of  the 
inverse  and  implicit  function  theorems  enables  us  to  obtain  the  main  result 
regarding  existence  of  approximate  solutions  and  error  estimates. 

Theorem  7.2:  (1)  Let  yQ  e  F(E),  X  =  Z©A  a  parameter  splitting  of  X, 

X  =  V  ©  T  a  coordinate  splitting  for  M  at  a  given  point  xQ  e  M,  and 


A:  Y  ■*  V  an  isomorphism.  Suppose  that  an  operator  (7.2)  and  a  basic  dis¬ 
cretization  are  chosen  so  that  the  stability  condition  (S)  is  satisfied. 
Then  for  all  sufficiently  small  h  >  0,  there  exist  xQh  e  such  that 

lim  x  .  =  x  . 
h->0  oh  0 

(2)  In  addition,  assume  that  r  >_  2  and  let  x:  JCT+M  be  a 
Cr-function  representing  the  manifold  M  locally  near  xQ  as  given  by 
Theorem  2.2.  Then  there  exist  a  compact  ball  ^cj,  0  e  JQ,  and  Cr- 
functions  x^.  JQ  representing  the  approximate  manifolds  ^  locally 

near  x0^  such  that 

(7.9)  | |x(t)-xh(t) | |  <  C||(I-Ph)Qx(t)||,  t  c  J0, 

where  C  is  independent  of  h  and  t. 


Actually,  we  can  say  more.  It  turns  out  that,  for  sufficiently  small  h, 
X  =  V©T  is  also  a  coordinate  splitting  of  X  for  each  at  x  .. 

Moreover,  if  x:  J  -*■  M  is  written  as  x(t)  =  xQ  +  t  +  An(t),  then  each 

xh:  Jo  +  ^  has  the  ^orm  xh^^  =  x0  +  t  +  Anh(t).  The  proof  of  Theorem  7.2 

follows  verbatim  the  corresponding  proof  in  [4  ]. 


8.  Discrete  Convergence  and  Stability  Conditions 


A  central  requirement  of  the  theory  in  the  previous  section  is  the 
condition  (S).  We  formulated  (S)  as  a  stability  condition  similar  to  those 
arising  typically  in  convergence  studies  of  discretization  methods.  In 
certain  practical  situations,  (S)  is  readily  verified.  For  instance,  in  [4] 
we  used  decompositions  F(x)  =  N(x)  +  G(x)  with  Q  =  DN(xq)  and  compact 
DG(xq),  in  which  case  the  stability  condition  (S)  holds  for  any  coordinate 
splitting. 

Nevertheless,  the  need  remains  for  other  necessary  and  sufficient  conditions 
for  the  validity  of  (S)  and  for  a  closer  analysis  of  the  interplay  between 
the  quantities  F,  Ph,  Q,  and  A  entering  into  (S).  This  is  the  topic  of  the 
present  section. 

We  follow  Vainikko  Cl 33  in  the  following  formulations  of  various  types 
of  discrete  convergence.  A  sequence  <xn>  of  elements  in  a  Banach  space  X 


is  called  discretely  compact 


for  short)  if  every  subsequence  of 


<xn>  has  a  convergent  subsequence. 

Let  (h  >  0)  and  B  be  bounded  linear  operators  from  a  Banach  space 
X  to  a  Banach  space  Y  such  that  B^x  -*•  Bx  for  all  x  e  X. 

The  convergence  B^  B  is  regular,  denoted  by  B^  £  B,  if 

(8.1)  <xh  :  hn  -*•  0>  bounded,  <Bh  xh  >  d-compact  =>  <x^  >  d-compact. 
"n  n  '  n  "n 

The  convergence  Bh  -*■  B  is  stable,  denoted  by  B,  ^  B,  if 


(8.2)  B’  exist  and  are  uniformly  bounded  for  sufficiently  small  h. 
h 

The  convergence  Bh  -*■  B  is  compact,  denoted  by  Bh  $  B,  if 


(8.3) 


<Xl  :  h  -*■  0>  bounded  *>  <Bh  xh  >  d-compact. 


Looking  back  at  our  condition  (S),  we  note  that 


0Ph(xo}Ay  =  (I-Ph)QAy  +  PhDF(x0)Ay  -  DF(x0)Ay,  y  e  Y, 
and  that  (S)  is  equivalent  with  stable  convergence  -  that  is,  with 
(8.4)  DFh(xQ)A  *  DF(xo)A* 

A  less  obvious  equivalence  is  the  content  of  the  following  result. 

Proposition  8.1 :  The  convergence  condition  (8.4)  and  hence  the  condition  (S) 
holds  if  and  only  if  Dfh(*0)A  £  DF(xQ)A. 

The  proof  follows  directly  from  a  theorem  of  Vainikko  [  1 3;  p.  655]  and  the 
observation  that  rge  DF(xQ)A  =  Y,  ker  DF(xQ)A  =  {0},  and  each  DFh ( xq) A 
is  a  Fredholm  operator  of  index  0. 

Both  regular  and  stable  convergence  involve  a  complex  interplay  among 
F,  P^»  Q,  and  A.  The  following  sufficient  condition  for  (S)  is  helpful 
here. 

Proposition  8.2:  If  Ph(-Q+DF(x0))A  $■  (-Q+DF(xq))A,  then  DFh(xQ)A  ^  DF(xq)A 
and  hence  (S)  holds. 

Proof:  From  DFh(xQ)A  =  QA  +  Ph(-Q+DF(xQ))A,  it  follows  that  QA  is  a 
Fredholm  operator  of  index  0.  If  QA  is  an  isomorphism,  then  we  have  trivially 
QA  £  QA  and  rge  QA  =  Y,  and  the  result  follows  immediately  [13;  p.  654].  If 
QA  Is  not  an  Isomorphism,  then  by  the  Fredholm  nature  of  QA  and  with 
Yq  ■  ker  QA,  Y^  *  rge  QA  there  exist  closed  subspaces  Y^  and  Y^  of  Y 
such  that  Y  *  Y  ©  Y,  *  Y2  ©  Y,,  dim  Y  *  dim  Y«  <  *  and  QA  is  an 


isomorphism  of  Yj  onto  Y y  Set  VQ  =  AYq  =  VfiA  and  V1  =  AY1 ;  then 
Vs*  VQ  ©  V^.  Let  =  X  ©  (Vj  ©  A)  and  9  =  VQ  ©  V^.  Since  dim  V  =  dim 
there  exists  an  isomorphism  B  of  VQ  onto  7  and  with  the  projection  P 
of  Y  onto  Yq  along  Y^  we  may  define  the  isomorphism 

A:  Y  ->  V,  Ay  =  A(I-P)y  +  BAPy,  y  e  Y. 

Then  QA  is  an  isomorphism  of  Y  and 

DFh(xQ)A  =  QA  +  Q(A-BA)P  +  Ph(-Q+DF(xQ))A. 

We  now  have  trivially  QA  £  QA  and  rge  QA  =  Y,  and  thus 

Q(A-BA)  +  Ph(-Q+DF(xQ))A  $  Q(A-BA)P  +  (-Q+DF(xq))A. 

The  result  follows  again  as  above  [13;  p.  654]. 

As  suggested,  the  converse  of  Proposition  8.2-is  not  true.  For  example 
let  Y  be  an  infinite-dimensional  separable  Hilbert  space  and  with  X  =  Yx 
consider  the  function 

F:  X  -*>  Y,  F(x)  =  |  y,  x  =  (y,A)  e  X. 

For  a  coordinate  splitting  X  =  V  ©  T,  take  V  =  Y  x  {0}  and  T  =  {0}  x  ]R, 
and  define  A:  Y  -*■  V  by  Ay  =  (y »0)  for  y  e  Y.  The  operator  Q:  X  -*•  Y 
is  defined  as  Qx  =  y  for  x  =  (y , A )  e  X  and,  if  {e.. }  is  an  orthonormal 
basis  for  Y,  the  discretization  projections  P^  are  taken  to  be  the  orthog 
onal  projections  onto  the  subspaces  spanned  by  e^....,e  with  —y  <  h  < 
Then  we  have,  for  any  x  , 


but 


Dph(xo)A  =  1  +  ?  ph  5  7  1  =  DF<xo)A’ 

Ph(-Q+DF(X0))A  -  \  ph  f\  1  =  (-Q+DF(xq))A. 

In  applications,  the  choice  of  Q  is  at  our  disposal  and  a  suitable 
choice  of  Q  may  ensure  the  validity  of  (S).  In  our  example,  the  slightly 

3 

different  operator  Q  defined  by  Qx  =  j  Y  f°r  *  =  (y.M  e  x  turns  out 
to  allow  the  application  of  Proposition  8.2. 

Compact  convergence  does  not  appear  to  be  any  easier  to  verify  than 
either  regular  or  stable  convergence,  and  so  Proposition  8.2  mainly  has 
theoretical  interest.  Our  next  result  is  a  step  towards  rectifying  that. 

Proposition  8.3: 

(i)  If  (-Q+DF(xq))A  is  a  compact  operator,  then 
Ph(-Q+DF(xo))A  $  (-Q+DF(x0))A. 

( 1 i )  If  Ph(-Q+DF(xQ))A  $  (-Q+DF(xq))A  and  if  Y  is  separable, 
then  (-Q+DF(xq))A  is  a  compact  operator. 

Proof:  Set  C  *  (-Q+DF(xq) )A.  To  prove  (i),  let  <yh  >  be  a  bounded 

n 

sequence  in  Y.  Since  C  is  compact,  any  subsequence  <Cyh  >  of  <Cyh  > 

nk  n 

has  a  convergent  subsequence  <Cy.  >.  If  y  is  the  limit,  then 
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where  |  i !  I  £a  for  sufficiently  small  h  (by  the  uniform  boundedness 

principle).  Consequently,  any  subsequence  of  <Ph  Cyh  >  has  a  convergent 

n  n 

subsequence  and  hence  P^C  *  C.  (ii)  follows  directly  from  a  result  of 
Vainikko  [13;  p.  654]. 


Proposition  8. 3( i i )  remains  true  for  certain  nonseparable  spaces  Y  as  well. 
However,  most  applications  involve  separable  spaces  and  so  there  is  no  real 
loss  in  stating  the  result  as  we  did. 

Proposition  8.3  effectively  replaces  the  need  for  testing  the  discrete 
convergence  by  a  simple  compactness  condition  on  the  operator  (-Q+DF(xo))A 
and  thus  represents  a  significant  simplification.  The  projections  P^  are 
removed  from  consideration,  and  the  condition  appears  to  involve  only  A  and 
Q.  We  may  actually  go  one  step  further  and  remove  A  from  the  picture.  The 
final  condition,  expressed  in  the  next  proposition,  is  a  condition  on  Q  and 
Q  alone. 


Proposition  8.4:  ( -Q+DF(xq ) ) A  is  compact  if  and  only  if  -Q+DF(xq)  is 

compact. 


Proof:  Set  K  *  -Q+DF(xq).  If  K  is  compact,  obviously  KA  is  compact. 

So  assume  KA  is  compact.  Let  <xn>  be  a  bounded  sequence  in  X.  Since 

X  =  V©T,  we  can  write  x„  =  v„  +  t  and  the  sequences  <v  >  and  <t  >  of 

n  n  n  n  n 


components  are  bounded.  Since  dim  T  =  m,  <tp>  has  a  convergent  subsequence 


>. 


Set  yn  =  A"  vn;  then  <yn  > 

k 


is  a  bounded  sequence  and  hence,  by  the 


compactness  of  KA,  <KAyri  >  has  a  convergent  subsequence  <KAy_  >.  Now  it 

nk  nk. 

follows  that 


so  that  <Kxn  >  converges.  Therefore,  K  is  compact. 
ki 

We  summarize  all  the  above  implications  and  equivalences  in  a  theorem. 

Theorem  8.5:  Let  yQ  e  F(E),  X  =  Z  €>  A  a  parameter  splitting  of  X, 

X  =  V©T  a  coordinate  splitting  for  M  at  a  given  point  xQ  e  M,  and 
A:  Y  -*■  V  an  isomorphism.  Suppose  that  an  operator  Q  (7.2)  and  a  basic 
discretization  are  chosen  and  consider  the  following  statements: 

(i)  stability  condition  (S)  is  satisfied; 

(ii)  DFh(xQ)A  ^  DF(xq)A; 

(iii)  DFh(xQ)A  5  DF(xo)A; 

(iv)  Ph(-Q+DF(x0))A  S  (-Q+DF(x0))A; 

(v)  (-Q+DF(xq))A  is  a  compact  operator; 

(vi)  -Q+DF(xq)  is  a  compact  operator. 

Then  the  following  implications  hold: 

(i)  <=>  (ii)  <=>  (iii)  <=  (iv)  <=  (v)  <=>  (vi). 

If,  in  addition,  Y  is  separable,  (iv)  <=>  (v). 

Thus,  condition  (vi)  becomes  an  extremely  simple  condition  on  the 
operator  Q  alone  which  guarantees  the  stability  condition  (S).  Since  (iii) 
and  (iv)  are  not  equivalent,  (vi)  is  not  the  only  way  to  the  stability 
condition  (S),  but  it  is  one  of  the  easiest  to  apply. 
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9.  Applications 

In  this  final  section,  we  look  at  two  applications.  The  first  example 
is  a  mildly  nonlinear  problem,  whereas  the  second  one  is  much  more  general 
in  nature. 


a.  Shallow  Arch  Subjected  to  a  Load:  This  problem  concerns  the  deformations 
of  a  shallow  circular  elastic  arch  which  has  been  used  frequently  as  a  test 
case,  and  we  refer  specifically  to  [7].  In  dimensionless  form,  the 
total  potential  energy  of  the  system  is  given  by 


(9.1) 


V  = 


{[v*  -  u  +  j-a^u')2]2  +  ct2(u")2  -  2a3p}  de. 


-0_ 


where  u  and  v  are  (dimensionless)  displacements  of  the  arch  axis,  the 
primes  denote  differentiation  with  respect  to  the  angle  0,  20q  is  the  angle 
subtended  by  the  circular  arph,  p  is  the  (dimensionless)  load,  and  aj,a2,a3 
are  (dimensionless)  constants.  The  load  p  =  p(0,u,v,X)  is  permitted  to  be 
a  sufficiently  differentiate  function  of  the  variables  0,u,v  and  of  m 
parameters  X^,...,Xm  which  we  represent  as  an  m-vector  X.  For  clamped 
ends,  the  boundary  conditions  are 


(9.2)  u(+0Q)  =  u'(±0o)  =  v(+0Q)  =  0. 

The  principle  of  stationary  potential  energy  applied  to  (9.1)  yields  the 
pair  of  differential  equations 

a2u,,n  -  a1a3Pvu'  -  [v'  -  u  +  ^  (u'  J2]^ u”  +  1)  -  a3Pu  =  0, 


v"  -  u'  +  a, u'u"  +  a  o,  =  0. 
1  3  * 


(9.3) 


These  equations  (9.3)  together  with  the  boundary  conditions  (9.2)  may  be 
formulated  weakly  as  the  problem  of  finding  (u,v)  e  H^C-6q »0o3  *  [-eQ , 0Q ] 
such  that 

<u,<*)>2  +  <f(u,v,A),4»o  =  0,  4>  e  H^[-eo>eoL 

(9.4) 

+  <g(u,v,x),^>0  =  o,  ip  e  hJ[-0o10o], 

where  f,g:  H2[-0O,0O]  x  hJ[-0o,0q]  x  lRm  -*■  L2[-0o,0q]  are  given  by 

f(u,v,A)  =  -  {a1a3Pyu'  +  [v1  -  u  +  j  (u* )2](a1u"  +  1)  +  c^}, 

g(u,v,A)  =  u'  -  o^u'u"  -  a3Pv, 
and 

0Q 

«p,\p>.  =  V1^*  i  =  0,1 ,2. 

1 

-eo 

Suppose  that  p  is  such  that  f  and  g  are  Cr-mappings  with  r  >  2. 

-2  2 

We  introduce  linear  operators  k:  H  [-0.0  ]  H  [-0  ,0  ]  and 

00  0  00 

l:  H_1[-0O,0O]  ->  [“©Q » ©Q-]  defined  by 

<ku ,4>>2  -  u  c  H  [“OqjQqI,  <t>  c  ^q[”®q>®q3> 

<£v ,!/<>•,  =  <v,ip>0,  v  £  H'1[-0O,0O],  ip  £  hJ[-0o,0q]. 

Then  k:  L2[-0O,0Q]  -*•  Hot"e0’e0^  and  l:  ^  Ho^‘eo,eo^  are  comPact 

and  (9.4)  becomes  simply 


(9.5) 


u  +  kf(u,v,A)  =  0, 


v  +  £g(u,v,A)  =  0. 


n  -i 

If  we  now  set  Y  =  H*[-eo,eo]  *  %[-e  ,e  ],  X  =  Y  x  IRm,  and  define 
K:  H“2[-0o,0o]  x  H'1[-0o,0o]  -+  Y  by  K(u,v)  =  (ku ,iv)  and 
G:  X  -*■  L2[-0q,0o]  x  L2[-0Q,eo]  by  G(u,v,A)  =  (f (u,v,A) ,g(u,v,A) ) ,  then 

(9.5)  reduces  further  to 


(9.6)  F(w,A)  e  w  +  KG(w,X)  =0,  w  =  (u,v)  e  Y,  X  e  IRm. 


Thus,  the  shallow-arch  problem  has  been  transformed  into  a  mildly  non¬ 
linear  problem  (9.6),  and  a  mildly  nonlinear  problem  can  be  handled  in  a 
straightforward  manner.  With  X  and  Y  as  above,  there  is  a  natural 
parameter  splitting  X  =  Z  ©  A,  where  Z  =  Y  x  {0}  and  A  =  {0}  x  ]Rm. 
Consequently,  a  natural  choice  for  the  operator  Q  e  L(X,Y)  of  (7.2)  is  the 
projection 


Q(w,X)  =  w,  x  ~  (w,A)  e  X. 


In  terms  of  Q,  (9.6)  can  be  written  as 


F(x)  =  Qx  +  KG(x)  =0,  x  e  X, 


and  hence  -Q  +  DF(xQ)  =  KDG(x0)  is  compact  for  any  xQ.  By  Theorem  8.5, 
the  stability  condition  (S)  is  satisfied  for  any  coordinate  splitting  X  =  V 
and  isomophism  A:  Y  V,  and  the  estimate  (7.9)  of  Theorem  7.2  becomes 


|w(t)-wh(t)MY  +  |  |X(t)-Xh(t)  |  |  jpm  <  C  ||(I-Ph)w(t)||Y,  t  e  JQ. 


Setting  P^w  =  (P^u.P^v),  we  may  write  this  estimate  in  terms  of  the 
original  displacements  u  and  v  as 

I |u(t)-uh(t) | j  2+  | |v(t)-vh(t) | |  1  +  | |X(t)-Xh(t) | |  jpm 
Ho  Ho 

1  c  {||(I-P|J1))u(t)l|  2  +  ||(I-P|S2))v(t)||  1),  t  e  J0. 

H  H 

0  o 

b.  A  Nonlinear  Pi  rich! et  Problem:  Let  ft  be  a  bounded  domain  in  IR^  with 
sufficiently  smooth  boundary.  As  usual,  we  use  the  multi-index  notation  Dau 
for  the  partial  derivatives  of  functions  u  defined  on  ft.  If  |a|  =  £a.  =  l 
and  there  is  no  danger  of  confusion,  we  write  also  l/u  for  the  generic 
derivatives  Dau  of  order  |a|  *  l. 

Let 

2k 

F(u,A)  5  $(C,u,Du,. . . ,D  u,A)  =  0 

be  a  given  elliptic  differential  operator  on  ft.  Here,  ?  e  is  the  space 
variable  and  A  e  IRm  a  parameter  vector,  and  $  denotes  a  sufficiently 
differentiable  function  of  all  its  arguments.  Given  any  appropriate  function 
f  defined  on  ft,  we  consider  the  nonlinear  Dirichlet  problem 

(9.7)  F(u,A)  =  f  on  ft , 

(9.8)  Dau  =0  on  9ft  for  0  £  |a|  <  k-1. 

Unlike  the  previous  problem,  which  was  transformed  into  a  mildly  nonlinear 
problem,  we  analyze  (9.7)  and  (9.8)  directly.  For  a  fixed  scalar  x,  0  <  x  <  1, 
let  W  be  the  Banach  space  of  functions  in  C  ’  (ft)  which  satisfy  the 
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boundary  conditions  (9.8).  Then  set  X  =  W  x  ]Rm  and  Y  =  C°*X(Q).  With 
this,  the  Dirichlet  problem  (9.7)-(9.8)  becomes  simply 

(9.9)  F(x)  =  f,  x  =  (u,X)  e  X. 

Again,  we  have  a  natural  parameter  splitting  X  =  Z  ®  A  with  Z  =  W  *  {0} 
and  A  =  {0}  x  Df. 

For  xQ  =  ( uQ , XQ ) ,  the  derivative  DF(xQ)  is  given  by 


Now  assume  that  xQ  is  a  solution  of  (9.9)  where  F(uq,Xo)  is  strongly 
elliptic;  that  is,  where  the  partial  derivative  DuF(xQ)  is  a  strongly 
elliptic  linear  differential  operator.  The  principal  part 


(9.10) 


Lx  =  ,  I  ~-d~  (g.un»Du  ,...,D2ku  ,X  )Da,  x=  (u,X)  e  X, 
Xo  |a|=2k  3(Dau)  00  00 


of  DuF(xQ)  is  then  an  isomorphism  of  W  onto  Y.  Moreover,  the  remaining 
terms  comprising  DuF(xQ)  are  all  compact  operators.  For  a  proof  of  these 


facts,  we  refer  to  [  1  ;  p.  686  ff.].  Hence,  it  follows  that  DuF(xq)  is  a 


Fredholm  operator  from  W  to  Y  with  index  0.  Let  Pw 
projections  of  X  onto  W  and  IRm,  respectively.  Then 


and 


be  the 


Pl(  is  a  Fredholm 
W 


operator  of  index  m  and  in  the  decomposition 


DF(x0)  -  DuF(x0)Pw  +  0AF(xo)PRm 

the  term  F ( xQ ) P  m  is  compact.  This  shows  that  DF(xq)  is  a  Fredholm 

1R 

operator  from  X  to  Y  of  index  m.  Therefore,  if  we  assume  that  xQ  is 

a  regular  solution  of  (9.9),  our  theory  applies. 

The  above  discussion  already  suggests  a  choice  for  the  operator 

Q  e  L(X,Y)  of  (7.2),  namely,  Q  =  L  .  In  fact,  then  -Q  +  DF(x  )  is 

x° 

compact  and  by  Theorem  8.5  the  stability  condition  (S)  is  satisfied  for  any 
coordinate  splitting  X  =  V©T  and  isomorphism  A:  Y  ->•  V.  The  estimate 
(7.9)  of  Theorem  7.2  now  takes  the  form 


2k,T+  l|X(t>-Ah(t)M  <  c  ||(I-Ph)L  u(t)||  0>T  .  t 


where  L  is  the  principal  part  given  in  (9.10). 


41 


References 


[1]  S.  Agmon,  A.  Douglis  and  L.  Nirenberg,  Estimates  near  the  Boundary  for 
Solutions  of  Elliptic  Partial  Differential  Equations  Satisfying  General 
Boundary  Conditions  I,  Comm.  Pure  Appl .  Math.,  12,  1959,  623-727. 

[2]  Yu.  G.  Borisovich,  V.  G.  Zvyagin  and  Yu.  I.  Sapronov,  Nonlinear  Fredholm 
Maps  and  the  Leray-Schauder  Theory,  Russian  Math.  Surveys,  32:4,  1977, 

1-54. 

[3]  F.  Brezzi,  J.  Rappaz  and  P.-A.  Raviart,  Finite-Dimensional  Approximation 
of  Nonlinear  Problems. 

Part  I:  Branches  of  Nonsingular  Solutions,  Numer.  Math.,  36,  1980,  1-25. 

Part  II:  Limit  Points,  Numer.  Math.,  37,  1981,  1-28. 

Part  III:  Simple  Bifurcation  Points,  Numer.  Math.,  38,  1981,  1-30. 

[4]  J.  P.  Fink  and  W.  C.  Rheinboldt,  On  the  Discretization  Error  of  Parametrized 
Nonlinear  Equations,  Univ.  of  Pittsburgh,  Inst,  for  Comp.  Math,  and 
Applications,  Tech.  Report  ICMA-82-40,  July  1982,  to  appear  in  SIAM  J. 

Numer.  Anal. 

[5]  J.  P.  Fink  and  W.  C.  Rheinboldt,  On  the  Error  Behavior  of  the  Reduced-Basis 
Technique  for  Nonlinear  Finite- Element  Approximations,  Zeit.  f.  Angew. 

Math,  und  Mech. 

[6]  V.  Girault  and  P.-A.  Raviart,  An  Analysis  of  Upwind  Schemes  for  the  Navier- 
Stokes  Equations,  SIAM  J.  Numer.  Anal.,  19,  1982,  312-333. 

[7]  A.  D.  Kerr  and  M.  T.  Soifer,  The  Linearization  of  the  Prebuckling  State  and 
Its  Effect  on  the  Determined  Instability  Loads,  Trans.  ASME,  J.  Appl.  Mech., 
36,  1969,  775-783. 

[8]  S.  Lang,  Differential  Manifolds,  Addison- Wesley,  Reading,  Massachusetts, 
1972. 

[9]  T.  Poston  and  I.  Stewart,  Catastrophe  Theory  and  Its  Applications,  Pitman, 
London,  1978. 

[10]  M.  Schechter,  Principles  of  Functional  Analysis,  Academic  Press,  New  York, 
1971. 

[11]  M.  Spivak,  A  Comprehensive  Introduction  to  Differential  Geometry,  Vol.  I, 

2nd  ed..  Publish  or  Perish,  Berkeley,  California,  1979. 

[12]  M.  M.  Vainberg  and  V.  A.  Trenogin,  The  Methods  of  Lyapunov  and  Schmidt  in 
the  Theory  of  Nonlinear  Equations  and  Their  Further  Development,  Russian 
Math.  Surveys,  17,  1962,  1-60. 

[13]  G.  Vainikko,  Approximative  Methods  for  Nonlinear  Equations  (Two  Approaches 
to  the  Convergence  Problem),  Nonlinear  Analysis,  Theory,  Methods  & 
Applications,  2,  1978,  647-687. 


